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ABSTRACT 


The purpose of this thesis is to design a quantum well infrared detector with 
tunable spectral bandwidth. The tunability of the bandwidth is achieved by using the 
linear Stark effect for the ground to first excited state transition in an asymmetric 
quantum well. The position of the absorption peak is dependent on the direction of the 
electric field, and therefore it can be either blue or red shifted by changing the direction 
of the field. If two identical asymmetric quantum wells are arranged opposite each other, 
we can obtain both the blue and red shift for either direction of the bias. This method can 
produce broader peaks with tunable bandwidths proportional to the applied field. 

A program was developed to calculate the energy levels and wavefunctions of an 
arbitrary quantum well. The program was used to design a step quantum well capable of 
detecting infrared in the 8-12 p ,m band. The validity of the approach was verified by 
comparison with experimental data of Martinet et al. (Martinet, 1992), and found to have 
a good agreement. The designed step well was used to create a tunable bandwidth 
detector. The analysis showed that the bandwidth could be tuned to more than twice the 
peak width. The numerical simulation indicates the possibility of manufacturing a tunable 
bandwidth infrared detector by using step quantum wells. 
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I. INTRODUCTION 


A. INFRARED DETECTORS 

The human eye can detect only a small segment of the electromagnetic spectrum, 
the familiar visible band ranging from blue to red. Only objects hot enough to radiate at 
these frequencies and those that reflect visible light are visible to the eye. For the most 
part, however, objects are too cold to glow noticeably, making them undetectable at 
night. Yet their finite temperature does give them an infrared glow. It is this glow that 
night-vision infrared cameras strive to see. At the heart of an infrared camera is its “eye”: 
a two-dimensional array of detector pixels, with each pixel converting some of the 
incident photons to an electric signal. 

A Quantum Well Infrared Photodetector (QWIP) consists of a stack of quantum 
wells designed to detect infrared light (Levine, 1993). When the quantum well is 
adequately deep and narrow, its energy states are quantized. The potential depth and 
width of the well can be adjusted so that it has only two energy states: a ground state near 
the well bottom and a first excited state near the well top. In Figure 1.1, photoexcited 
electron transition from the ground state to the first excited state, where an externally 
applied voltage sweeps them out, producing a photocurrent. 



Figure 1.1 Electron excitation of a quantum well infrared photodetector. 

For two bound states in the quantum well, only photons having energies 
corresponding to the energy separation between the two states are absorbed, resulting in a 
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detector with a sharp absorption spectrum. If there is only one bound state in the well, the 
transition occurs between the bound state and the continuum, producing a broader peak 
(Levine, 1993). Designing a quantum well to detect light of a particular wavelength 
becomes a straightforward matter of adjusting the potential depth and width of the well to 
produce two states, separated by the desired photon energy. Impressive advances in the 
crystal growth and processing techniques of compound semiconductors such as gallium 
arsenide ( GaAs ) and aluminum gallium arsenide ( Al x Ga x _ x As ) with different aluminum 
compositions allowed the fabrication of novel QWIP detectors (Levine, 1993) and arrays 
(Gunapala, 1996). 

B. PROGRESS OF QWIPS 

In the 1970's, InSb and HgCdTe were used to cover the 3-5 and 8-12 micrometer 
wavelength ranges. However, HgCdTe has some shortcomings due to difficulty in its 
material preparations, uniformity and reproducibility in the fabrication process. Other 
difficulties include mechanical softness, sensitivity to elevated temperature, and low yield 
for large array fabrications, which make such detectors extremely expensive. In order to 
fabricate low-cost detectors and overcome the above difficulties, detection using 
superlattices was suggested by Esaki and Tsu (Esaki, 1969). Extensive investigations, 
both theoretical (Coon, 1984 and 1985) and experimental (Smith, 1983), were a subject 
of considerable interest and were carried out until the first observation of strong 
intersubband transitions in GaAs / AlGaAs quantum wells by West and Englash (West, 
1985). Since the problems associated with material preparation were avoided, this new 
development has received considerable attention. 

The first bound-to-bound quantum well infrared photodetector was demonstrated 
by Levine et al. (Levine, 1987). Because the photoelectrons had to tunnel out of wells 
from the bound excited state under large barriers, as illustrated in Figure 1.2(a), the 
responsivity was so low that the device had to work using a C0 2 laser. This low 
responsivity led to the use of bound-to-continuum transitions, which were first proposed 
by Coon and Karunasiri in 1984 (Coon, 1984). Hasnain et al. (Hasnain, 1989) 


2 



demonstrated the first bound-to-continuum QWIP (see Figure 1.2b) to have a much 
higher responsivity, but also a small peak absorbance due to broad spectral width. 



Figure 1.2 Schematic diagrams showing (a) bound-to-bound and (b) bound-to- 

continuum electron transfer. 

Several configurations of the QWIP structure have been reported based on 
transitions from bound-to-extended states (Levine, 1990), bound-to-quasi-continuum 
states (Levine, 1991), bound-to-quasibound states (Gunapala, 1996), and bound-to- 
miniband states (Yu, 1991). These configurations provide good detection performance 
because in all designs, electrons in the ground state do not flow in response to the bias, 
while electrons activated to the excited state yield photocurrent under low bias. Several 
structural optimizations have also been made to improve the performance of QWIPs such 
as dark current reduction from thermionic emission (Gunapala, 1996). 

C. PURPOSE OF THIS THESIS 

The goal of this research is to design a quantum well detector with a tunable 
spectral bandwidth. The tunability of the bandwidth is achieved by using “the linear Stark 
effect” in an asymmetric quantum well (Yuh, 1989). For example, in the step quantum 
well shown in Figure 1.3(a), the amount of Stark shift from the ground to first excited 

eFci 

states transition is approximately ——, where F is the electric field strength, and a is 
the step width (Yuh 1989). 
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Figure 1.3 Schematic diagrams showing an (a) asymmetric step quantum well and a 
(b) designed step quantum well structure for bandwidth tuning. 

As previously mentioned, the position of the absorption peak is sensitive to the 
direction of the electric field and can be either blue or red shifted by changing the 
direction of the field (Mii, 1988). If two wells are arranged as in Figure 1.3(b), both red 
and blue shifts can occur at the same time producing an absorption peak with a tunable 
bandwidth. 

In the following chapter, we will present a detailed analysis of energy levels in a 
step quantum well, including external electric field effects. In addition, the absorption 
coefficient under applied field will be calculated using a step approximation to the linear 
potential. Finally, in Chapter 3, we will present an optimized device structure to achieve 
tunable bandwidth detection. 
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II. THEORETICAL BACKGROUND 


A. INTRODUCTION 

In the following sections, we will develop the theoretical tools necessary to 
analyze quantum well energy levels as a function of applied bias for an asymmetric 
quantum well. The Schrodinger equation is solved numerically using the transfer matrix 
approach (Wang, 1989). The effect of the electric field on the energy levels will be 
incorporated by the “staircase” approximation of the linear potential. 


B. TRANSFER MATRIX METHOD 

In the case of an infinite square quantum well, the Schrodinger equation is given 

by 


h 2 d 2 'V, 


= ETC 


n = 1,2,3,... 


( 2 . 1 ) 


2 m e dz 

where zis the growth direction, n is the quantum number, fi is the reduced Planck’s 
constant, m* e is the effective mass of the electron, and 'P B and E n are the wavefunction 


and energy in the n ,h energy level, respectively. The wave function ? must also satisfy 
boundary conditions across a potential step since the probability current across a 
boundary should be continuous. 


In light of the above the conditions and the fact that the wave functions go to zero 
at the boundaries (z=0,L) and the electron stays within the well, we can obtain the 
energy and wavefunction of the quantized energy states as 

E n = n 2 ——T-j n = 1,2,3,... (2.2) 

2 m e L 

and 


V P„ (z) = A sin 


( m z ^ 


n =1,2,3,... 


(2.3) 
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where L is the width of the quantum well. Taking the origin at the left side of the well as 
shown in Figure 2.1, the wave function can then be normalized over the total distance 
- °o < z < +°° as follows: 

]'¥:(z)'¥ n (z)dz = l, (2.4) 



Figure 2.1 Infinite square well potential and energy levels. 

There are many techniques available to calculate the energy levels and the wave- 
functions of finite quantum well structures, each having advantages and disadvantages. 
The deciding factor on which method to use depends on, to some extent, the structure 
itself. For example, the variational approach yields simple results for the ground state but 
has difficulty finding the excited states. The solution of the Schrodinger equation for a 
finite square quantum well can be found by matching boundary conditions (Weisbuch, 
1991). However, when the number of interfaces is large, this method becomes 
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cumbersome (Choi, 1997). With the help of a computer, the transfer-matrix method is a 
multipurpose approach in obtaining the energy states and the corresponding wave 
functions of a structure with an arbitrary potential profile. In the following section, the 
transfer matrix method will be discussed. 

A heterostructure is made from layers of materials having different bandgaps. For 
example, in a typical square quantum well, a layer of GaAs is surrounded between the 
two thicker higher bandgap layers of Al x Ga^_ x As . The well width is controlled by 
adjusting the GaAs layer thickness while the potential depth is controlled by adjusting 
the Al composition in the barrier layers. Again the application of the boundary 
conditions leads to the quantized energy values E n and wave functions 'P n (z). 


In the transfer matrix method, a plane wave is assumed to be incident from the left 
of the structure as shown in Figure 2.2. In each layer the wave function can be written as 

a linear combination of the forward ( A n e lkz ) and backward ( B n e lkz ) traveling plane 

waves. A matrix can be set up using the boundary conditions at each interface of the 
quantum well structure to connect the wavefunctions between the adjacent regions 
(Hutchings, 1989). 


V(z) 


Incoming Wave 


777 , 


777 , 


777 , 


777 


n— 2 


m , 


! 777 , ! 

I n -1 | 


Figure 2.2 Arbitrary multi-layered potential well structure. 

For the arbitrary quantum well structure shown in Figure 2.2, the wave function 
for the n th region can be written as 
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X) f n (z) = A n e i ^+B n e- ik ‘ 


n = 1,2,., N 


( 2 . 6 ) 


where A n and B n are arbitrary constants, 


K = ^jf(E-VJ, ( 2 - 7 ) 

m* n is the electron effective mass in the n" 1 layer and the thickness of the n' h layer 
d„ = z„ +l -z„. Furthermore, E is the electron energy, which can either be higher or 

lower than the highest potential in the structure. The boundary conditions at the n' h 
interface give (Robinett, 1997) 

V„(4)=V„+i(z„) 


1 dy„(z n ) _ 1 d\\r n+l (z n ) 

m l dz n dz n 

Thus, the following matrix, M n , can be obtained to connect the coefficients in the n lh 
layer to that of the (n + \) th layer as follows: 


where 


and 


rA ,i 


\A 1 

n+ 1 

= M„ 

n 

A+1_ 


A. 


n = 1,2 ,...,n -1 


M =- 

" 2 


(1+y )e l ^ kn kn +^ Zn (l-y )e l ^ k,i+k n+i) Zn 

^_y ^ e d k n +k n+\) z n ^ + y ^ e ~d k n~ k n+] ) z n 


y _ K m n + 1 

' n ; * 

k n + i m n+l 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 


In Equation (2.10) z n is the spatial coordinate of the interface between n' h and (/? + 1)'' 


layers. 


This model can be extended to the multi-layered quantum well structure and by 
using the above relationship, the coefficients of the wavefunctions for the first and last 
layers are related by 



m„m„_ 1 ...m 2 m 1 


A " 

= M 

A ' 

_ 

m n 

772,2 

' A ' 

A . 


A. 


_ m 21 

™ n _ 

A. 


( 2 . 12 ) 
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If the two outer layers have higher potentials than the energy of the particle (i.e., if 
E < V ), then there will be bound energy states. In this case, we can set the coefficient 
B n = 0, because it would become infinite as z —» 00 . In the same manner, A x can be set 
to zero for the region z < 0 . Under these conditions, Equation 2.12 reduces to 

(2.13) 


(2.14) 


"V 


n h 1 

m \2 

'o' 

_ 0 


_ m 21 

m 22 _ 

A. 


In order to satisfy the Equation 2.13, we must have 

m^E) = 0 . 


A Matlab program was developed to plot m 22 (E) versus E and identify the points where 
m 22 is zero. For example, Figure 2.3 shows the calculated m 12 as a function of E for a 
square potential well. The calculated bound state energies and wavefunctions are shown 
in Figure 2.4. 


x -jo*’ Finding Energy Levels In a Potential Wei 



o 

Figure 2.3 Plot of m 22 (E) versus E for a square potential well with L = 110 A, 

V X =V 3 = 0.75 eF and V 2 = 0 eV . 
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Figure 2.4 Energy levels and wavefunctions of the quantum well described 

in Figure 2.3. 


C. INTERSUBBAND TRANSITIONS 

In a quantum well structure made of semiconductors, the wavefunction consists of 
both envelope and Bloch functions. Thus, not only the transitions between Bloch 
functions but also between the envelope states are allowed within the structure. The 
transitions between Bloch states correspond to interband transitions, while transitions 
between envelope states correspond to intersubband transitions. For a given interaction 
potential V , the transition rate W can be obtained from the Fermi’s golden rule (Coop 

1984) 

1T = — IF,,IT*, >l 2 5 (Ef-Ei-tm), (2.15) 

h f 
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where E i and E f are the electron energies in the initial and final states, respectively, and 
ft co is the incident photon energy. The interaction potential energy, V , in the dipole 
approximation is given by (Coon, 1984) 


C = 


in 


in 


n 1/2 


2e 0 n r (pc 


e- p 


(2.16) 


where e is the electron charge, I is the incident photon flux, c is the speed of light, £ 0 
is the electric permittivity, n r is the index of refraction, (0 is the angular frequency of the 


incident photon, e is the unit vector of the photon polarization direction, and p is the 
momentum operator, 


p = —ifi 


d d d ^ 

i -h j -I -k — 

v c)x dv dz J 


(2.17) 


In the case of a quantum well, the wavefunctions for the initial and final states of 
transition can be written as 


¥, = «,■( r)^,(z) 


(2.18) 


where n ( and u f are the Bloch functions of the initial and final states, respectively, 


k i and k f are the initial and final wave vectors in the xy plane, r is the 3-D position 


vector, p is the position vector in the xy plane, TV and 4' are the corresponding 

envelope functions, and A is the area of the structure. The transition matrix element in 
Equation 2.15 can be evaluated as (Coon, 1984) 
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<^ / \V P I'V, >=<u f (r^ f (z)^\V p \u i (r^ i (z)^> 


=<uAr)e' kf ' P IV \u i (r)e ik ‘ p x\|/ f (z)l\|/,(z)> (2.19) 


+ <y f ( z )^j =-1 V p I \| i. (z) e -jj >< u f (r) I m. (r) >. 


In the above expression, we assume that the envelope functions vary slowly compared to 

—> —> 

the Bloch functions. The integration over the x and y coordinates gives k i = k f . The first 
term of Equation 2.19 is zero because of the fact that the Hamiltonian used to obtain \|/ ( . 
and \\i f is Hermitian. In addition, the Bloch functions are normalized to unity. Thus, the 
inter subband transition matrix element reduces to (Coon, 1984) 


<'i’ f \V p \'i’ i >=<y f (z)\V p \y i (z)> 


JL 

in 


in 


n1/2 


( 2 . 20 ) 


v 2e 0 n r G) J 


<\l/ f (z)l/ 3 J\)/,(z)>cose, 


where 0 is the angle between the photon polarization angle and the z axis and P, is the 


momentum operator in the z direction. The strength of the transition is usually 
characterized by the oscillator strength, which is defined as 


where 


/ = 


2m* (E 2 ~E t ) 


h- 


(Vf Izl’*'/) 


( 2 . 21 ) 


('P / lzl'P I .) =j 'i’fZ'i’idz. 


( 2 . 22 ) 


D. INTERSUBBAND TRANSITION SELECTION RULES 

The first selection rule of intersubband transition is the parity selection rule. From 
Equation 2.20, it is found that the parity of the final state has to be different from that of 
the initial state in order to obtain a non- zero transition matrix element. This means that 
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the intersubband transitions in a square quantum well can occur only from odd states to 
even states or vice versa. For example 1 —» 2 or 2 —> 3 transitions are allowed whereas 
1—>3 is not allowed. Figure 2.5 shows schematically the parity selection rule of the 
intersubband transitions in a square potential well. However, the symmetry of a square 
well can be broken by applying a relatively strong electric field, making the normally 
“forbidden” transitions become allowed (Mii, 1990). 


barrier 


n 

X 


4 


4 


e 2 

E, 


0 


a 


Figure 2.5 Square quantum well showing allowed and forbidden transitions. 


The second selection rule is the polarization selection rule. Since the potential 
modulation of a quantum well structure is only in one direction, the z direction, 
intersubband transition will not occur when the photon polarization is perpendicular to 
the z direction. This can be understood by considering the cos0 term in Equation 2.20. 
This unique polarization selection rule is an important characteristic and can be used for 
the identification of the intersubband transition (West, 1985). In order to obtain a strong 
intersubband transition, the photon polarization has to coincide with the growth direction 
(i.e., 0 =0). 

For an asymmetric quantum well (see Figure 2.6), the dipole matrix elements of 
transitions for both odd-to-odd and even-to-even quantum numbers do not vanish, since 
the eigenfunctions of these energy states are not eigenfunctions of the parity operator. 
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Figure 2.6 Asymmetric quantum well with allowed transitions. 


In addition to allowing intersubband transitions between any low states, the 
asymmetric well structure also gives linear Stark shift (Yuh 1989). This can be utilized 
for tuning the detection wavelength (Martinet, 1992). 


E. INTERSUBBAND ABSORPTION COEFFICIENT 

The absorption coefficient in quantum well devices is one of the most important 
factor in the design of the quantum structures and is proportional to the amount of 
incident light that is absorbed in the structure. The absorption coefficient depends upon 
the material characteristics of the respective semiconductor. For example, the high 
refractive index of a semiconductor such as GaAs , (n = 3.3), causes that light incident on 
the surface of an epitaxial layer structure is refracted towards the normal, for example, 
and away from the required orientation. Efficient coupling of the intersubband transition 
can be achieved by a number of optical techniques, including the use of a beveled edged 
device (Levine, 1987), diffraction gratings (Hasnain, 1989), and combined diffraction 
gratings and waveguides (Andersson, 1992). 

In addition to absorption, it is important to extract the photoexcited carriers out of 
the quantum well. This is achieved by placing the excited state near the top of the barrier 
as shown in Figure 2.7. In the case of the bound-to-continuum transition (see Figure 
2.7a), the excited state is located just above the barriers, while for a bound-to-quasi 
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bound structure (see Figure 2.7b) the excited state is just below the barrier. The latter 
structure provides lower leakage current due to the higher barrier height for the same 
detection wavelength. 




(a) 




(b) 


Figure 2.7 Bound-to-continuum (a) and (b) bound-to-bound transitions. 

For the bound-to-bound transitions, since there is only one final state, the 
summation over the initial states in Equation 2.23 gives the number of the electrons in the 
ground state, 


N = N d LA . (2.23) 

Here, N is the number of electrons in the initial state, N d is the three dimensional 

doping density, and L and A are the well width and area of the quantum structure, 
respectively. The delta function in Equation 2.15 is usually replaced by the normalized 
Lorentzian function to take into account the spectral line width given by (Yariv, 1997) 


g(t> co) = 


j_r_ 

2* (E f -E-tmf + (y 2 ) 2 ' 


(2.24) 
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where T is the full width at half the maximum. This is usually obtained experimentally 
and is about 10-20 meV for intersubband absorption. The transition rate, found using the 
interaction potential and the momentum operator is given by 


' Bound —> Bound 


N d (IAL)e 2 ti 3 
(hCd)[m* e ) 8 0 n r c 


k V/ I —ly,>l 2 g(ftco). 


(2.25) 


For the bound-to-continuum transitions, the number of initial states is still N 
with energy E i , as in the case of bound-to-bound states. However, the summation over 
the final states requires multiplication by the joint density of states p (E), which is 


defined as 


p {E)dE- 


(2.26) 


For energy states separated by the incident photon energy, the above equation can be 


reduced to 


P(£) 27 Ah 2 (E f -V b ) 


(2.27) 


In this expression, / is the length of the region used for the normalization of the extended 
state wavefunction, m* b is the electron effective mass in the barrier, and V b is the 
potential height of the barrier. The transition rate for bound-to-continuum states is then 


defined as 


r Bound —» Continuum 


N d {IAL)e 1 ti 1 l 2m* d 2 

—/ - - - —^- \<w f \ — lur. >1“. 

2 (m,) £,n r c i(*«» (£f~K) ' * 


(2.28) 


This model of the bound to free intersubband absorption is one where the line width is 
due only to the variation of the joint density of states. The absorption coefficient a(ftco) 

can be obtained by using both the transition of the matrix element and the transition rate 
as follows: 


a (tun) = 


number of transitions per unit volume and time „ , /at 

- ; — ; ---, a 0oo) = x 


(2.29) 


incident photon flux I 

Thus, the absorption coefficients for the bound-to-bound and bound-to-continuum 
transitions are given by 
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2*3 


a 


Bound > Bound 


(7?co) = 


N d e n 


(TK£>){m* e ) e 0 n r c 


I<¥/ I-7-IV/ >|2 S( n ®) 
dz 


(2.30) 


a 


Bound —>Continuum 


(tm) 


X c hi 


2m, 


2(m e ) 8 0 n r c^®) (^7 


l<\|/, I — l\j/, >l 2 


f dz 


The length term, /, which appears in the bound-to-continuum transition, is eliminated by 
the term in the normalized final state wave function (assuming it is confined to a 


length, /). As an example, absorption coefficients for the bound-to-bound and bound-to- 
continuum situations for the quantum well shown in Figure 2.8 are depicted in Figures 
2.9 and 2.10, respectively. 


Quantum Wall with Zaro External Electric Field 
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Figure 2.8 Energy levels of the Al 040 Ga 060 As / Al 020 Ga 0m As /GaAs / Al 040 Ga 0m As 

quantum potential well. 
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Figure 2.9 Bound-to-bound absorption coefficient for the quantum potential well 

given in Figure 2.8 with T - 20 meV . 


Figure 2.9 shows the absorption versus photon energy for bound-to-bound 
transitions. There is a sharp rise in the absorption spectrum when the photon energy is 
equal to the difference between the bound state energies in the quantum well structure. 
The width of the curve depends on the finite lifetime of the excited electron, which was 
modeled using the Gaussian spectral density function. 
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Bound to Continuum Absorption 



Figure 2.10 Bound-to-continuum absorption coefficient for the quantum well 

in Figure 2.8. 

Similarly, Figure 2.10 shows the absorption coefficient versus the photon energy 
for the case of bound-to-continuum transitions. The sharp rise in the spectrum 
corresponds to a photon energy equal to the difference between the bound state energy 
and the energy at the top of the quantum well where the extended continuum states begin. 
As the photon energy increases above this threshold value, absorption increases until 
reaching a peak value corresponding to the energy above the top of the barrier. At very 
high energies there is a slow decrease of a (fib)) since the density of states drops as the 

energy is increased. 

As previously mentioned, our goal is to design a detector with tunable bandwidth. 
The most convenient way to do this is by using bias as the tuning parameter. As a result, 
it is important to study the effect of bias on the energy levels as well as on the absorption. 
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F. EFFECT OF ELECTRIC FIELD ON QUANTUM WELLS 

When a uniform electric field is applied along the z direction (growth direction), 
the quantum well structure will tilt in the opposite direction of the electric field and the 
change due to this tilt will affect the energy states of the quantum well. The potential 
experienced by an electron in the quantum well is given by 

V(z)=V n +eFz, (2.31) 

where V n is the potential height of the n' 1 ' layer of the quantum well, e is the electron’s 
charge, F , measured in (V / m ), is the electric field strength , which is applied along the 
growth direction. A square quantum well in the presence of an external bias is shown in 
Figure 2.11. 



Figure 2.11 Schematic representation of a square quantum well under an applied bias. 

By increasing the electric field, the lowest energy state E x in the quantum well 

will decrease in energy and fall into to the triangular potential formed at the bottom of the 
well. On the other hand, the excited states are less affected by the external electric field. 
However, for low electric field strengths, the shifting of all the energy levels are the same 
due to the symmetry of the square well. Thus, there is no appreciable Stark effect. This 
can be altered using an asymmetric quantum well, which is equivalent to a square well 
under a high electric field. 
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1. Transfer Matrix Approach for Quantum Wells Under Applied Bias 

Previously, we have discussed the calculation of the energy levels and 
wavefunctions of arbitrary quantum wells using transfer matrix approach. However, this 
approach cannot be used directly due to the linear variation of the potential energy under 
an external bias. 

One way to overcome this problem is to approximate the linear potential due to 
the bias shown in Figure 2.12, with a series of steps or “staircase” as shown in Figure 
2.13. The advantage of this approach compared to the use of Airy functions is that, it can 
be used at low field strengths, whereas Airy functions tend to diverge at low bias 
voltages, making the numerical analysis much more difficult (Brennan, 1987). 


Quantum Wei with 5x10 6 V/m External Electric FieW 



Figure 2.12 Step quantum well with applied bias of 5 V /|im . 
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Figure 2.13 The stair case approximation of the potential shown in Figure 2.12. 

The Schrodinger equation for an electron in the quantum well in the presence of 
an electric field F can be written as 

“T-FtT 'V,-(E-V,-qFz)'V,=0, (2.32) 

2 m n dz 

where E is the energy of the electron in the layer and V n and m n are the potential and the 


effective mass in the n"' layer, respectively. In order to apply the transfer matrix relations 
it is necessary to know the boundary conditions at the end points. For bound states, we 
can use the same approach as used in the zero bias case. For the continuum states we can 
assume that the electron is incident from the left side of the potential shown in Figure 


2.13, and which gives B N = 0. Thus, starting from 



all the unknown 


1X1 


1X1 


coefficients can be progressively calculated, ending with the coefficient 

N 

_ B n_ 

0 


Therefore the electron wave functions can easily be obtained because all the coefficients 
are produced when \ is known (Hutchings, 1989). 


22 



Figures 2.14 (a) through 2.14 (d) show the potentials and eigenenergies calculated 
using the “staircase” approximation for external bias varied from 0 to 5 VI \\m. In this 
way, we can observe the change in the energy states as a result of the applied electric 
field. Table 2.1 shows the change in bound state energies due to the external applied bias 
in the quantum well. 




Quartn M «t*xl Eitami Eftrtic Rec 


-4 4 0 


E i =2* T it*V 


»J3«n - 


< • I 10 11 M (0 

«mi ■«* 


(a) 


aurtr rtM nrVIT r>*™ Brctc hm 



QmviM* 





Ej*2SmV 














' t>’5n 

-4-2 0 7 

4 « 1 M 

m 

(b) 

12 « it 

•n' 

OMjn VM m St Wi Eon* ftett AM 



in* 


■ N' 


(c) 


(d) 


Figure 2.14 Effect of electric field on the bound state energies of a step quantum well 
for the following field strengths: (a) F = 0, (b) F = 1 VI [ xtn , (c) F = 3 VI (im and 

(d) F = 5 Vl[im. 
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Zero 

IE/ |im 

3 V / (im 

5 V / (im 

Energy States 

Biased Well 

Biased Well 

Biased Well 

Biased Well 


(meV) 

(meV) 

(meV) 

(meV) 



47 

61 

74 

Ej 

41 

?E=6 

?E=20 

?E=33 



153 

164 

175 

e 2 

147 

?E=6 

?E=17 

?E=28 



250 

256 

— 

e 3 

247 

?E=3 

?E=9 

Continuum 


Table 2.1 Bound state energies as a function of bias. 


2. Intersubband Absorption Under an Applied Electric Field 


The intersubband absorption peak position will shift in energy with changing 
electric field due to the shift of the energy levels. At a low bias, this effect is especially 
pronounced in asymmetric quantum wells as compared to symmetrical wells. Figure 2.15 
shows the calculated absorption coefficient as a function of bias across the well using 
Equation 2.30. The matrix element in Equation 2.30 is evaluated using the wavefunction 
under the applied bias. 
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3000 


Bound to Bound Absorption With Applied Bias 



Figure 2.15 Absorption coefficient as a function of bias across the step quantum well 

structure. 

The strength of absorption is found to increase with increasing electric field, 
mainly due to the stronger matrix element. The asymmetry observed between positive 
and negative biases is due to the fact that under forward bias (negative), the first excited 
energy level moves toward the ground state and hence the transition energy is reduced, 
while under reverse bias (positive) the opposite is true and the transition energy is 
increased. 

G. COMPARISON WITH EXPERIMENT 

In order to verify the validity of the “staircase” model for a quantum well under 
bias, we have analyzed a step quantum well structure employed by Martinet et al. 
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(Martinet, 1992) for the study of bias dependent photocurrent. The parameters of the 
structure used in Martinet’s experiment are shown in the Table 2.2. 


Entered Experimental Data 

Layer 

and Thickness 

Concentration 

AlGaAs (barrier) 
36.5 run 

0.40 

AlGaAs (step) 6.7 

0.20 

GaAs 

3.2 run 

— 

AlGaAs (barrier) 
36.5 nm 

0.40 


Table 2.2 Parameters of the step quantum well structure used by Martinet et al. 

(Martinet, 1992). 

Figure 2.16 shows the step well structure and the calculated energy levels using 
the Matlab program developed using the “staircase” approximation for the same structure 
used by Martinet et al. (Martinet, 1992). 
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Figure 2.16 Calculated energy levels of the Al 040 Ga Q60 As / Al 020 Ga 0m As IGaAs step 

quantum well. 


The calculated value of E 2 — E l of 107 meV is close to the zero-field absorption 
peak of 115 meV observed experimentally. The experimentally observed shift of the 
absorption peak (Martinet, 1992) for different electric fields in the step quantum well is 
shown in Figure 2.17. 
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Photoresponse (normalized) 



Energy (meV) 

Figure 2.17 Normalized responsitivity of the asymmetric step quantum well under 
different positive and negative biases. The peak shift is nearly four times the spectral 
width. Zero field absorption (bold line) is included as a reference, (Martinet, 1992). 


Normalized Program Output for Different Electric Fields 



Figure 2.18 Plot of bound-to-bound, 1—>2 normalized absorption transitions with 

r = 20 me V done with the Matlab program. 
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Figure 2.18 shows the calculated absorption coefficient (normalized) for the same 
set of bias voltages as in the experiment. The experiment shows that the peak position 
changes from 93 meV to 145 meV when the bias is changed from 42 kV /cm to 
-40 kV /cm . The calculated value of the peak changes from 86 meV to 125 meV for 
the same bias values. The experimental and calculated peak shifts amount to 52 meV and 
39 meV, respectively. The difference between the two values may be due to the 

comparison of the photocurrent versus the absorption coefficient. Nonetheless, the results 
show that the model predictions are in reasonably good agreement with the experimental 
observation, indicating the accuracy of the calculation. 

H. SUMMARY 

In this chapter the transfer matrix method was derived for an arbitrary potential 
well. This technique was employed to calculate the energy levels and wavefunctions of a 
step potential well. Intersubband transition selection rules in both symmetric and 
asymmetric potential wells and how they affect the design of the detector were also 
discussed. Another significant design parameter, the absorption coefficient, was 
calculated for both symmetric and asymmetric wells. External bias and its effects on the 
quantum wells, energy levels, and absorption coefficients were discussed. The use of 
“staircase” approach to calculate the absorption coefficient is found to be in good 
agreement with the experimental observation. In the following chapter, the design of the 
tunable bandwidth quantum well detector will be discussed. 
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HI. DESIGN OF A TUNABLE BANDWIDTH AlGaAs/GaAs STEP 
QUANTUM WELL INFRARED PHOTODETECTOR 


A. INTRODUCTION 

As discussed in Chapter II, intersubband absorption results from transitions 
between the energy levels of a quantum well formed within the conduction or valence 
band. This section presents detailed calculations of energy level, wavefunction and 
oscillator strength in an AlGaAs / GaAs step well structure. The use of an electric field to 
tune the detection bandwidth will also be described. First, however, the calculations of 
effective mass and barrier height as a function of Al composition in the AljGa x _ x As layer 
will be discussed. 

B. EFFECTIVE MASS AND BARRIER HEIGHT CALCULATION 

Normally the band offsets of the conduction band (A E c ) and valence band ( AE ) 
of Al x Ga l _ x As / GaAs quantum well structures can be calculated by using the total band 
offset (A E ). Accepted values of conduction and valence band offset ratios are as 
follows (Meyer, 2001): 

AE = 0.6 AE 

(3.1) 

AE v = 0.4 AE g . 

The bandgap of Al x Ga x _ x As , E g (x) , can be approximately calculated using Vegard’s 
Law as (Adachi, 1994) 

E g (x) = xE g(AlAs) + (1 — x)E g (Ga4i ) , (3.2) 

where x is the Al molar fraction and E g(AIAs) and E g(GaAs) are the bandgaps of AlAs and 
GaAs , respectively, at the T - point. 

Using the experimentally observed band gaps of AlAs (2.671 eV) and 
GaAs (1.424 eV) , the band gap of Al x Ga t _ x As at the T - point is given by (Singh, 1993) 

E g (x) = 1.424+1.247 x eV (3.3) 
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The graph in Figure 3.1 shows the band gap of the Al x Ga^ x As layer as a function of 
Al composition where the composition x changes from zero to one hundred percent. 


Energy Gap vs. Molar Fraction of Al 



Figure 3.1 Energy gap at T - point as a function of Al. 

1. Effective Mass Calculation 


The effective mass of the Al x Ga^ x As alloy can also be found using Vegard’s 

Law as 


1 

* 

m (AIGaAs ) 


X 


m (AlAs) 


+ 


\ — x 

m (GaAs) 


(3.4) 


The inverse relationship is due to the fact that the effective mass and energy are inversely 
related, where 


m (AlAs) 


= 0.15 m e 


n \GoAs ) =0.067 m e 


(3.5) 
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are the T point effective masses of AlAs and GaAs, respectively(Singh, 1993). 
Figure 3.2 shows the plot of effective mass as a function of Al composition. 


Effective Mass vs Molar Fraction of Al 



Figure 3.2 Effective mass vs. molar fraction of A/. 


33 



THIS PAGE INTENTIONALLY LEFT BLANK 


34 



IV. FINAL DESIGN OF THE STRUCTURE 


In the following, we will discuss the design of a tunable bandwidth detector with 
center wavelength near 8.5 \im . The most convenient material system is AlGaAs/ GaAs 
due to their lattice match and to mature growth technology. In order to demonstrate the 
band tuning capability, a bound-to-bound transition in a step quantum well will be used. 
Figure 4.1 shows a step well with peak absorption near 8.9 [im . This structure was 
optimized to achieve a relatively large absorption coefficient with peak absorption 
between 8 and 9 \lm . The first excited state lies about 24 meV below the barrier edge. 


Quantum Well without External Electric Field 



Figure 4.1 Optimized quantum well showing the energy levels, centroids, and 

lambda. 

Figure 4.2 shows the calculated absorption coefficient of the structure for a set of positive 
and negative bias voltages. 
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Bound to Bound Absorption With Applied Bias 



Figure 4.2 Absorption coefficient for bound-to-bound intersubband transition for the 
step well shown in Figure 4.1 for a set of bias voltages. 


Table 4.1 shows E x , E 2 , and oscillator strength, / , for the quantum well structure given 
in Figure 4.1 under various applied biases. 


Energy states 
and oscillator 
strength 

Zero 

Biased Well 

IF/ \Xm 

2 V / \im 

3 V t ]Xm 

V 

-V 

V 

-V 

V 

-V 

E,(eV) 

0.0766 

0.0833 

0.0767 

0.0926 

0.0792 

0.1017 

0.0818 

E 2 (eV) 

0.2160 

0.2184 

0.2220 

0.2211 

0.2275 

0.2224 

0.2334 

f 

0.1343 

0.1146 

0.1523 

0.0960 

0.1831 

0.0766 

0.1798 


Table 4.1 Energy states and oscillator strengths under different biases for the 
quantum well structure shown in Figure 4.1. 
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It can be seen in Figure 4.2 that the absorption peak shifts with applied bias due to 
the linear Stark effect in the step quantum well. However, the bandwidth of the 
absorption peak remains nearly the same. In order to achieve tunable bandwidth, a two 
step quantum well was arranged as shown in Figure 4.3. 


Quantum Weil 'without External Electric Field 



Figure 4.3 Final design of the asymmetric quantum well structure. 


In this case, under a given bias the absorption peaks of the two wells move in 
opposite directions. Figure 4.4 shows the calculated absorption coefficient of the 
combined structure as a function of photon energy for a set of bias voltages. 
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Bound to Bound Absorption With Applied Bras 



Figure 4.4 Calculated absorption coefficients for the combined step quantum well. 


It can clearly be seen that the bandwidth increases as the bias increases. At 10 
kV / cm , the bandwidth is nearly doubled. For higher electric fields, the absorption 
coefficient near the center drops due to a large shift of the two peaks. The details of the 
layer structure needed to fabricate the tunable bandwidth quantum well detector are 
shown in Table 4.2. 
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Layer 

and Thickness 

Molar Fraction 

Layer 

m* (kg) 

E g (eV) 

AlGaAs (barrier) 

5 nm 

0.32 

Alo. 32 Gao. 68 As 

0.0814mo 

1.8230 

AlGaAs (step) 

7 nm 

0.25 

Alo. 25 Gao. 75 As 

0.0778mo 

1.7357 

GaAs 

4.8 nm 

-- 

GaAs 

0.067mo 

1.4240 

AlGaAs (barrier) 

5 nm 

0.32 

Alo. 32 Gao. 68 As 

0.0814mo 

1.8230 

GaAs 

4.8 nm 

-- 

GaAs 

0.067mo 

1.4240 

AlGaAs (step) 

7 nm 

0.25 

Alo. 25 Gao. 75 As 

0.0778mo 

1.7357 

AlGaAs (barrier) 

5 nm 

0.32 

Alo. 32 Gao. 68 As 

0.0814mo 

1.8230 

Table 4.2 The layered structure of a tunable bandwid 

1 th step quantum well detector. 


In summary, the design of a step quantum well structure for tunable bandwidth 
detection was introduced. The tunability of the detector by applying a bias to the quantum 
structure was also demonstrated. 
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V. CONCLUSION 


The transfer matrix method was derived, and energy levels and wavefunctions for 
an asymmetric arbitrary potential well were calculated. Intersubband transition selection 
rules and their effect on the design process were also discussed. The absorption 
coefficient, which is also important in the design of a quantum well, was calculated for 
both symmetric and asymmetric wells. The “staircase” potential method was used to 
present both the effects of an external bias and the shift in the detection wavelength due 
to the various applied external electric fields across the detector. External bias and its 
effects on the quantum wells, energy levels, and absorption coefficients were explained in 
detail. The validity of the program was demonstrated through comparison with 
experimental data and the best possible design of the step quantum well structure for 
broadband detection was found. 

A quantum well structure having a broadband detection capability was designed 
using two step quantum well structures arranged opposite to each other. This kind of 
structure allowed both blue and red Stark shifts of the transition energy producing a 
tunable bandwidth. The calculated absorption spectral bandwidth as a function of the 
applied bias shows that the peak width can be varied from 1.27 p m to 2.17 pm , which 
is twice the bandwidth at zero bias. 

Finally, the analysis indicates the possibility of fabricating a tunable bandwidth 
infrared detector using asymmetric quantum wells. 
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APPENDIX. [SELECTED PROGRAMS] 


%THIS PROGRAM FINDS BOUND ENERGY LEVELS IN ANY KIND OF STEP WELL 
%by ATAKAN KONUKBAY © 2002 
clear all,clc,close all; 
global N x El h Ein v 

N=input('Enter the number layers> '); % ENTER 1 FOR SQUARE,2 FOR STEP QUANTUM WELL 

e=1.602e-19; %Coulomb 
h=197.33e-9; % Planck's Constant 

%MOLAR RATIOS (ENTER VALUES AS .35) 

XX(l)=input('Enter the percentage of A1 in region I >'); 

XX(2)=input('Enter the percentage of A1 in region II >'); 

XX(3)=input('Enter the percentage of A1 in region IV >'); 

%THE ENERGY GAP IN EACH REGION OF THE ASYMMETRIC QUANTUM WELL 

Eg 1=1.424+1.247*XX( 1); % AlxGa( 1- x) As 

Eg2=l.424+1.247*XX(2); %AlxGa(l-x)As 

Eg3=1.424; %GaAs 

Eg4= 1.424+1.247 *XX( 1); % AlxGa( 1- x) As 

%EFFECTIVE MASSES 
M=.511e6; 

m(l)=l/((XX(l)/(0.15*M))+((l-XX(l))/(0.067*M))); 

m(2)=l/((XX(2)/(0.15*M))+((l-XX(2))/(0.067*M))); 

m(3)=0.067*M; 

m(4)=l/((XX(3)/(0.15*M))+((l-XX(3))/(0.067*M))); 

% POTENTIAL IN CONDUCTION BAND 
v(l)=0.6*(Eg 1-Eg3); 

v(2)=0.6*(Eg2-Eg3); %Set origin @ 0 
v(3)=0; %.6*(Eg3-Eg2); 

v(4)=0.6*(Eg 1-Eg3); 

cy c * *** ** 


a=input('Enter the width from 0->a in Angstrom> '); 
b=input('Enter the width from a->b in Angstrom> '); 

x=[0 a a+b]*le-10; 

cy c * *** ** 


E=0.0001; 

Ein=input('Enter the energy (eV)> '); 

Efin=input('Enter the external E field (le6 V/m)> '); %ENTER AS 1,2...ETC. 
El_=Ein*le4; 
while E<=Ein; 
for d=l:El_ 

T=l; 

for n = 1 :N+1 

k(n)=sqrt((2*m(n))*(E-v(n))/((h A 2))); 

k(n+l)=sqrt((2*m(n+l))*(E-v(n+l))/((h A 2))); 
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A(n)=(exp(i*(x(n)*(k(n)-k(n+l)))))*((k(n+l)*(m(n)))+(k(n))*(m(n+l)))/(2*(k(n+l))*(m(n))); 
B(n)=(exp(-(i*(x(n)*(k(n+l)+k(n))))))*((k(n+l)*(m(n))) -(k(n))*(m(n+l)))/(2*(k(n+l))*(m(n))); 
C(n)=(exp(i*(x(n)*(k(n)+k(n+l)))))*((k(n+l)*(m(n)))-(k(n))*(m(n+f)))/(2*(k(n+f))*(m(n))); 
D(n)=(exp(-(i*(x(n)*(k(n)-k(n+l))))))*((k(n+l)*(m(n)))+(k(n))*(m(n+l)))/(2*(k(n+l))*(m(n))); 


MM(n)={[A(n) B(n);C(n) D(n)]}; 
end 


ZZ=1; %Unit matrix 
for n=l:N+l; 

ZZ=MM{n}*ZZ; 

end 

M22=real(ZZ(2,2)); 

TT(d)= {[T*M22,E]}; 

Tt(d,:)=TT{d}; 

E=E+.0001; 

end 

end 


zz=l; 

for z=l :E1_-1 
tti(z)=Tt(z,i); 
tt2(z+l)=Tt(z+l,f); 
ttm=ttf(z)*tt2(z+l); 
if ttm<0 
E(zz)=[Tt(z,2)]; 
zz=zz+l; 

end 

end 

end 

clc 

MLDRAWFINE 

%MLDRAWFINE 

%THIS PROGRAM DRAWS THE SHAPE OF THE QUANTUM WELL & DIVIDES THE WIDTH OF 
THE STRUCTURE INTO MANY STEPS 
%by ATAKAN KONUKBAY © 2002 

clc; 

e= 1.602e-19;%Coulomb 
Edp=Efin; %Bias 
delta=le-10; 

X{ 1,1 }=(x(l)-50e-10):delta:x(l); 
for n=2:N+l 

X{ l,n}=x(n-f):delta:x(n); 
end 

X{ 1 ,n+1} =x( 1 ,n):delta: (x( 1 ,n)+50e -10); 
for n=l:N+2 
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XL(l,n)=length(X{ l,n}); 

Y{l,n}=(Edp*X{ l,n})*le6;%cnst; 
YL( 1 ,n)=length( Y {1 ,n}); 
end 

for n=l:N+2 

lgth=length( Y {n}); 

Y_{ l,n}=[Y{l,n}(l:lgth)]+v(n); 

end 


x_=[]; 

y—=[]; 

for n=l:N+2 

x_=union(x_,X {n}); 
y_=cat(2,y_,Y_{n}); 

end 

for n=f :N+2 
hold on; 

plot(X{ l,n},(Y {l,n}+v(l,n))); 
end 


for n=2;N+2 
LXf=XL(f,n); 

LYf=YL(f,n); 
hold on 

line([X{n-l }(l,XL(l,n-l)),X{n-l}(l,XL(l,n-l))],[Y{n-l}(l,YL(l,n-l))+v(n-l),Y{n}(l,l)+v(n)]) 
hold on 
end 

for n=l:length(E) 
switch N 
case f 
hold on 

line([x(l),x(2)],[E(n),E(n)],'Color','r') 
hold on 
case 2 
if E(n)<v(2) 
hold on 

line([x(2),x(3)],[E(n),E(n)],'Color','r') 

hold on 

else 

line([x(l),x(3)],[E(n),E(n)],'Color','r') 

end 

end 

end 

title(sprintf('Quantum Well With External Electric Field %g 10 A 6 V/m',Efin)); 

for n=l:N+2 
QA{ n }=Y {n} +v(n); 

L(n)=length(QA{n}); 

end 
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switch N 


case 1 

AQ=QA{2}; 

AQ(1)=[]; 

AQ((L(2>1))=[]; 

QA{2 }=AQ; 
v_=[]; 
case 2 

AQ=QA{2}; 

AQ(i)=G; 

QA{2}=AQ; 

AQ=QA{3}; 

AQ(i)=G; 

QA{3}=AQ; 

v_=[]; 

end 

for n=l:N+2; 

v_=cat(2,v_,QA{n}); 

end 

L=length(v_); 

vl=v(l); 

xlabel(sprintf('a = %g nm , b = %g nm',x(2)*le9,(x(3)-x(2))*le9)); 
ylabel(sprintf('Energy (e V)' ,Efin)); 
set(findobj(gca,'Type','line','Color',[0 0 1]),... 

'Color','blue',... 

'LineWidth',1) 

text(x(3)+28e-10,0+.01 ,sprintf('%g / %g / %g',XX(l),XX(2),XX(3)),... 

'HorizontalAlignment','left','FontSize',8) 
for n=l:length(E) 

text(x(3)+le-10,E(n),sprintf('E%g=%g ev',n,E(n)),... 

’HorizontalAlignment','left','FontSize',8) 

end 

grid off; 

%THIS PROGRAM CALCULATES THE ENERGY LEVELS FOR MULTI STEP WELL 
%by ATAKAN KONUKBAY © 2002 

global N v x El h Ein 
N=length(x_)-1; 


e=1.602e-19; %Coulomb 
h=197.33e-9; %Planck’s constant 

x=x_; 

%****** m ass 

m_l(l:XL(l))=m(l); 
m_2( 1: (XL(2) -1 ))=m(2); 
m_3(l:(XL(3)-l))=m(3); 
m_4(l :(XL(4)))=m(4); 
mm=[m_l m_2 m_3 m_4]; 
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v=v_; 


E=0.0001; 

Edrw=0; 

Ein=input('Enter the energy level (eV)> '); 

Efin=0 
El_=Ein*le4; 
while E<=Ein; 
for d=l:El_ 

T=l; 

for n = 1:N+1 

k(n)=sqrt((2*mm(n))*(E-v(n))/((h A 2))); 
k(n+1 )=sqrt((2*mm(n+1 ))* (E-v(n+1 ))/((h A 2))); 

A(n)=(exp(i*(x(n)*(k(n)-k(n+l)))))*((k(n+l)*(mm(n)))+(k(n))*(mm(n+l)))/(2*(k(n+l))*(mm(n))); 

B(n)=(exp(-(i*(x(n)*(k(n+l)+k(n))))))*((k(n+l)*(mm(n)))- 

(k(n))*(mm(n+l)))/(2*(k(n+l))*(mm(n))); 

C(n)=(exp(i*(x(n)*(k(n)+k(n+l)))))*((k(n+l)*(mm(n)))-(k(n))*(mm(n+l)))/(2*(k(n+l))*(mm(n))); 

D(n)=(exp(-(i*(x(n)*(k(n)- 

k(n+l))))))*((k(n+l)*(mm(n)))+(k(n))*(mm(n+l)))/(2*(k(n+l))*(mm(n))); 

MM(n)={[A(n) B(n);C(n) D(n)]}; 
end 

ZZ=1; %Unit matrix 
for n=l:N+l; 

ZZ=MM{n} *ZZ; 
end 

M22=real(ZZ(2,2)); 

TT(d)={ [T*M22,E]}; 

Tt(d,:)=TT{d}; 

E=E+.0001; 

end 

end 


zz=l; 


for z=l:El_-l 
ttl(z)=Tt(z,f); 
tt2(z+ l)=Tt(z+ f, 1); 
ttm=tt f (z)*tt2(z+1); 
if ttm<0 
E(zz)=[Tt(z,2)]; 
Edrw(zz)=E(zz)+.001; 

zz=zz+l; 

end 

end 
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%THIS PROGRAM CALCULATES THE ALPHA FOR BOUND TO CONTINUUM TRANSITIONS 
%by ATAKAN KONUKBAY © 2002 


clc,close 

e=1.602e-19; %C 

Nd=lel8/le-6; %l/m3 

hbar=1.055e-34; %J.s 

EoE=8.85e - 12*sqrt( 13); %F/m 
c=3e8 %m/s 


ENG=v(l)+.0001:.0001:v(l)+.02;%-(E(2)-E(l)):.001:2*((E(2)-E(l))); 

Alpha=Nd*e A 2*hbar A 2*sqrt((2*m( 1)*9.1 le-31/.5 1 Ie6)./((((ENG-E( l)).*e). A 2).*((ENG- 
v(l)).*e)))./(2*((m(3)*9.11e-31/.511e6) A 2)*EoE*c)*AvXRegBC. A 2 

plot(ENG,Alpha/le2) 
grid on 

xlabel('h\prime\omega [eV]') 
ylabel('Absorption [cm A { -1}]') 
titleCBound to Continuum Absorption') 

%THIS PROGRAM CALCULATES THE ALPHA FOR BOUND TO BOUND TRANSITIONS 
%by ATAKAN KONUKBAY © 2002 

clc,close 

Nd=lel8/le-6 %l/m3 
e=1.602e-19; %C 
hbar=1.055e -34 %J.s 
T=20e-3*e/2 %J 
EoE=8.85e-12*sqrt(13); %F/m 
c=3e8 %m/s 

HW=E(l)+.0001: .0001 :E(2); 

G=T./((((E(2)-E(l)-HW)*e). A 2)+((T) A 2)); 

CB_B=(Nd*(e A 2)*(hbar A 3))/(((m(3)*9.1 le -31/.51 le6) A 2)*EoE*c); 

Alpha=(CB_B.*G.*AvXRegMB_B A 2)./(HW.*e) 

plot(HW, Alpha/1 e2) 
grid on 

xlabel('h\prime\omega [eV]') 
ylabel('Absorption [cm A { -1}]') 
titleCBound to Bound Absorption') 

%THIS PROGRAM CALCULATES THE WAVEFUNCTION FOR MULTISTEP WELL 
%by ATAKAN KONUKBAY © 2002 

clc; 

figure; 

delta=le-10; 


trans=MM; 
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trs=MM; 
for n=2:N+l 

trans{ 1 ,n} =trs{l,n}*trans{ l,n-l}; 
end 

A(1,1)=0; 

B(l,l)=l; 
for n=l:N+l 
wavtrans=trans {1 ,n}; 

A( 1 ,n+1 )=wavtrans( f ,2) * B (f, f); 

B(l,n+l)=wavtrans(2,2)*B(l,l); 

end 

A(l,f)=0; 

B(i,f)=i; 

B(f ,(XL(f)+XL(2)+XL(3))-l:n+l)=0; 

X{ f ,f }=-x(2):delta:x(f); 
for n=2:N+f 

X{ l,n}=x(n-l):delta:x(n); 
end 

X{ l,n+l }=x(l,n):delta:(x(l,n)+(x(l,n)-x(l,n-l))); 

Totalfntegral=0; 
for n=l:N+2 

W{ f,n}=A(f,n)*exp(i*k(f,n)*X{ f,n})+B(f,n)*exp(-(i*k(f,n)*X{ f,n})); 

%%%***for the wave used in calculating the oscillation with-B parameter 
W_BC{ 1 ,n}=A(l,n)*exp(i*k(l,n)*X{ l,n})-B(l,n)*exp(-(i*k(l,n)*X{ l,n})); 

%%%*** 


W_{ l,n}=conj(W{ l,n}); 

integr{ l,n}=(W{ l,n}.*W_{ l,n})*delta; 

integregnj l,n}=sum(integr{ 1 ,n}); 

Totallntegral=integregn{ l,n }+Totallntegral; 
end 

bl = l/(sqrt(TotalIntegral)); 

normwave=0; 

for n=l:N+2; 

WavRealf l,n}=bf *W{ l,n}; 

%%%***Now multiplying with normalization factor to get the wave 
WavRealBC {1 ,n} =b 1 * W_BC {1 ,n}; 

%%%*** 


WR_{ l,n}=conj(WavReal{ l,n}); 

lnteg{ l,n}=(WavReal{ l,n}.*WR_{ l,n})*delta; 

lntegral{ l,n}=sum(lnteg{ l,n}); 

normwave=lntegral {1 ,n} +normwave; 

end 

for n=l:N+2 
hold on; 

plot(X{ l,n},WavReal{ l,n}); 

hold on; 

end 
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